Magnetic field dependence of the penetration depth of d-wave 
superconductors with strong isolated impurities 

Shan- Wen Tsai* P. J. Hirschfeld 

Department of Physics, University of Florida, Gainesville FL 32611 



Abstract 

A d-wave superconductor with isolated strong non-magnetic impurities should exhibit an upturn in the penetration 
depth at low temperatures [1]. Here we calculate how an external magnetic field supresses this effect. 
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In a recent paper [1] we presented a calcula- 
tion of current response of a d-wave superconduc- 
tor containing a single impurity and showed that it 
is singular in the low-temperature limit, leading in 
the case of strong scattering to a 1 /T term in the 
penetration depth. For a small number of such im- 
purities, this low-T upturn could be observable in 
cuprate superconductors. An estimate of the size 
of this effect and of the temperature range in which 
it should occur agrees with experiments reported 
by Bonn et al. [2] in which upturns correlated with 
disorder introduced by Zn atoms in YBCO have 
been observed. Here we consider the influence of a 
magnetic field on this effect. 

Since the low-temperature upturn in the pene- 
tration depth is due to the large number of quasi- 
particle excitations near the nodal directions, as 
is the similar upturn in the case of Andreev sur- 
face states [3-6] , we might expect any physical ef- 
fect which smears the gap nodes to cut off the up- 
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turn. In particular, the orbital coupling to an ap- 
plied magnetic field (nonlinear electrodynamics) 
will suppress the upturn as it does in the Andreev 
case. We need to add the Doppler shift of the quasi- 
particle energy iio — > ito + v s • k and follow the 
same steps as in Ref. [1] to calculate the magnetic 
field dependence of the upturn. Here v s is the lo- 
cal superfluid velocity and a typical shift v s kpj_ ss 
(H/H n )A , where H is the applied field and H n = 
3$o/(t'" 2 CoAo) is of the order of the thermodynamic 
critical field. 

The Nambu propagator for a pure <i-wave su- 
perconductor in the presence of an external mag- 
netic field is G k (w„) = [(iuj n + v s • k)r + £ k r 3 + 
A k Ti]/D k , where L» k = (ui n - iw s • k) 2 + g + A 2 ., 
£ k = £k — Mi an d the Tj are Pauli matrices. 
In the presence of one single 5-function im- 
purity, the Green's function can be expressed 
exactly in terms of the Green's function of 
the pure system and the T-matrix for the im- 
purity. The To-component of the T-matrix is 
TqH = (7rA )- 1 Go(^)/(c 2 - G (lo) 2 ) where 



Preprint submitted to Elsevier Science 



2 February 2008 



> 



5, 

o 



\W A = ° 


ri =0.01 T 

'p c 


0.0 r\ 




0.03\ 




0.05 











0.01 



0.02 



0.03 



T/T 



Fig. 1. Normalized change in the penetration depth due to 
single isolated impurity vs. normalized temperature T/T c 
for U! 3 /A = 0, 0.01, 0.03, 0.05, 0.09. 

c is cotangent of the s-wave scattering phase 
shift (c = corresponds to infinitely strong 
scattering) and G = (l/27riV )Tr £ k G£(w). 
Here we consider the case in which the su- 
pcrffuid velocity is in the x-direction, so the 
Doppler shift v s • kp at the four nodes is wi ; 4 = 
— <^2,3 = uj s = v s kp/V2- The integrated Green's 



function G (u>) = ( 



)K(- 



A 



)K(- 



— ^ where w± = w ± iu> s 



and K is the complete elliptic integral of the first 
kind [7]. We now use Eq. 3 of Ref. [1] to calculate 
the paramagnetic response 5K xx (p, q) and from 
it obtain the change in the penetration depth 
6 X following the procedure as in Ref. [1]. Fig. 1 
shows the result for a unitary scatterer (c = 0) 
and lu s /A = 0,0.01,0.03,0.05,0.09. The integral 
over momentum is performed for p, q ~ 1/A <C kp 
and nodal expansion of the order parameter. The 
Matsubara sum has been performed numerically. 

In order for the upturn to be observable, a fi- 
nite number of isolated strong scatterers must be 
present in the sample. An estimate of the likeli- 
hood of occurence of these rare configurations for 
randomly distributed impurities in order to extend 
the theory to finite disorder can be found in Ref. 
[1] . In real systems clustering may take place and 
we expect dependence not only on the concentra- 
tion of impurities but also on the nature of the sta- 
tistical distribution of the impurities in the sam- 
ple. Indeed measurements on YBCO single crystals 



doped with the same amount (0.31%) of Zn show 
strong sample-to-sample dependence [2] . 

In the absence of the external magnetic field, the 
parameter r\ p — vp/X sets a lower bound for the 
temperature in which the upturn occurs. A rough 
estimate for the value of the external magnetic field 
in which supression of the upturn should occur is 
given by us s — r/ p , i.e., H / H ~ £,oz/X. Typical 
values for YBCO (£ 0z « bA, X « 15001, H w 
2.5T) give H w 8mT. 

Observation of this effect would represent a 
remarkable example of the strong influence that 
statistically rare impurity distributions can some- 
times have on the macroscopic properties of a 
condensed matter system. It can be distinguished 
from the similar effect in a ci-wave superconduc- 
tor due to Andreev surface states by its disorder 
dependence and by its independence of field ori- 
entation and crystal shape. 
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